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ABSTRACT 

Adopting the advection dominated accretion flow (ADAF) model, we have re-examined the 
radiation emitted by isolated stellar-mass black holes accreting from the interstellar medium 
of our Galaxy. Two distinct types of black holes are given consideration, 1) remnants of 
conventional massive star evolution (RBHs, mass ~ IOMq), and 2) low-mass objects which 
possibly constitutes the dark halo and cause the observed microlensing events (MACHOBHs, 
mass ^ 0.5Mq). In contrast to previous studies incorporating spherical accretion models, the 
characteristic hard X-ray emission from ADAFs serves as an important observational signature. 
Standard estimates for the present density of RBHs suggest that a fair number of objects, either 
in nearby molecular clouds (such as Orion) or in the general Galactic disk, may be detected 
by near future X-ray instruments. Their observation may also be feasible at lower frequencies, 
particularly in the infrared and optical bands. On the other hand, the smaller mass and larger 
velocity on average of MACHOBHs compared to RBHs make them unlikely to be observable in 
the foreseeable future. 



Subject headings: accretion, accretion disks — black hole physics — dark matter — infrared: 
general — ISM: general — X-rays: general 



1. Introduction 

Massive stars are believed to end their lives undergoing gravitational collapse, with the more massive 
ones leaving behind black holes. The existence of such stellar-mass black holes in Galactic X-ray binary 
systems, radiating through gas accretion from their companion stars, has achieved strong observational 
support (e.g. Tanaka and Lewin 1996). Stellar-mass black holes should also occur in isolation (remnant 
black holes, hereafter RBHs), although in the absence of a significant accretion fuel reservoir, their 
confirmation is expected to be a difficult task. However, for those RBHs lying in high density environments 
such as molecular clouds and possessing small velocities relative to the surrounding gas, the rate of accretion 
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directly from the ambient interstellar medium (ISM) may be sufficiently high for the resulting emission to 
be observable (Shvartsman 1971). 

More spcciilativcly, black holes may also originate from processes other than conventional stellar 
evolution. The discovery of microlensing events towards the Large Magellanic Cloud by the MACHO 
collaboration imply that a sizable fraction (0.62^0 J) of the dark matter in our Galactic halo may consist 
of compact objects (MACHOs) in the mass range O.S^q ^Mq (Alcock ct al. 1997). At present there is no 
consensus on the most likely identity of MACHOs, and a wide variety of proposals have been put forth. 
While none are firmly excluded, the simplest sort of baryonic MACHOs, such as low-mass stars or collapsed 
objects resulting from ordinary stellar evolution all face various theoretical and observational challenges 
(e.g. Carr 1994). A serious possibility may then be that MACHOs arc low-mass black holes (MACHOBHs) 
that were somehow created in the early universe; indeed, a few such scenarios have already been suggested 
(Yokoyama 1997, Jedamzik 1997). 

Our objective here is to examine if accretion from interstellar gas onto RBHs or MACHOBHs in our 
Galaxy may give rise to observable radiation and aid in their identification. Pregalactic black holes may 
also accrete from the surrounding medium and radiate at early epochs, but Carr (1979) has shown that for 
the MACHOBH mass range of our interest (mass ^ IMq), no important constraints are imposed from the 
integrated background light, even if a high radiative efiiciency is posited. Thus, our consideration will be 
restricted to black holes in the solar neighborhood. 

Most previous investigations of ISM-accreting black holes have employed spherical accretion models, 
such as that of Ipser & Price (1982, 1983) (hereafter IP). However, there is a multitude of alternative 
spherical accretion models in the literature (for a review, see e.g. Chakrabarti 1996), and the emission 
spectra are sensitive to assumptions concerning the physical state of the accreting plasma and the radiation 
processes involved. Since none of them have been subject to observational tests, the validity of the 
particular IP model is in question. The advent of the advection dominated accretion flow (ADAF) model 
(Ichimaru 1977, Rees et al. 1982, Narayan & Yi 1994, Abramowicz et al. 1995; see Narayan 1996, 1997 
for reviews), being relevant for systems with relatively low accretion rates, has created a recent surge 
of interest in low-luminosity black holes. With minimal assumptions and few free parameters, ADAFs 
self-consistently predict a stable, hot, two temperature structure that generates broad-band spectra from 
radio to gamma-rays, and have been successfully applied to a variety of low-luminosity objects. The 
conditions for isolated, ISM-accreting black holes may also be appropriate for ADAFs, and in view of the 
fair amount of observational support, wc adopt ADAF as the model of choice. Note that Mahadevan & 
Quataert (1997) have also commented on the implications of ADAFs for isolated black holes with ISM 
accretion. 



2. Emission from RBHs 

We first quote some standard results for accretion from a uniform medium. (Bondi & Hoyle 1944, 
Frank, King & Raine 1992). A black hole of mass M moving with velocity V with respect to the ambient 
medium captures gas within the accretion radius defined by 

2GM 



where G is the gravitational constant and Cg is the medium's sound speed. In the following we ignore Cg, 
which is justified so long as the gas is cold (as in molecular clouds) and does not lie in a hot phase of the 
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ISM. The accretion rate is 

^ U%/ V 102 cm-ay 1^10 km s-iy 
» 5^3 X l„-.A/„. (^) (^) (^^) (2) 

where pgas and rigas are respectively the mass and number density of the medium, and 
MEdd = 1.39 X 10^^ (M/Mq) g s-i is the Eddington accretion rate with unit efficiency. 

Whether the gas then accretes onto the black hole more or less spherically or through a disk-like 
configuration depends on the amount of angular momentum available in the ISM. One possible source 
is ISM turbulence, for which the turbulent velocity i^turb has been measured to correlate with scale r as 
^'turb ~ 1-1 (r/lpc)'' '^^km s~^, albeit with large observational errors (Larson 1981, Falgarone & Phillips 
1990). Adopting this relation, the criterion that a rotating accretion flow can extend out to radius 
with specific angular momentum /; times the Keplerian value is that the black hole velocity satisfies 
F ^ 52 (rg /rofi'^)°'^^ (M /MQ f-^^km s^'^ , where Vg = 2GM/c^ is the Schwarzschild radius (Ipser & Price 
1977). Inhomogeneties in the medium (which is ubiquitous in the ISM) may also lead to accretion with 
angular momentum (Ipser & Price 1977), although there are conflicting opinions on this matter (Davies 
& Pringle 1980, Anzer, Borner & Monaghan 1987). ADAFs actually allow self-consistent solutions even 
for sub-Keplerian outer boundary conditions (see below), and we believe that they represent the most 
legitimate mode of accretion from the ISM. 

The principal assumptions of the ADAF paradigm are: 1) a viscosity prescription (Shakura & Sunyaev 
1973); 2) viscous heat primarily channeled into protons; 3) proton-electron energy transfer mediated by 
Coulomb collisions; and 4) existence of magnetic fields near equipartition. Models have been constructed for 
black holes of IOMq, typical of RBHs, accreting at different rates (corresponding to different combinations 
of rigas and V , Eq. [g|) in the ISM. We exploit global solutions of the dynamical equations, which include the 
important effect of energy advection and compressional heating for electrons (Nakamura et al. 1996, 1997a, 
Esin et al. 1997). The prescription for obtaining the solutions and computing the emission spectra are 
identical to that described in Manmoto, Kusunose & Mineshige (1997), except here the angular momentum 
has been set to half the Keplerian value (/; = 0.5) at To = lOOOr^ Q The free parameters in the model 
are chosen to be a = 0.1 for the dimensionless viscosity, and (3 = 0.5 for the ratio of gas pressure to total 
pressure. 

The results are shown in Figure 1. In contrast to the IP spherical accretion model or the standard 
thin disk model of Shakura & Sunyaev (1973), the spectrum arising from the hot, optically thin ADAF is 
very broad band and consists of roughly three peaks. Going from low to high frequencies, the processes 
responsible for each are synchrotron radiation (radio to optical), single Compton upscattering (UV to soft 
X-ray), and bremsstrahlung augmented by multiple Compton upscattering (hard X-ray to soft gamma-ray). 
Since hard X-rays are little affected by interstellar absorption, the occurrence of the high-energy component 
is a crucial point for the observability of RBHs. 



^Strictly speaking, this outer boundary condition conforms to the above ISM turbulence relation only for objects with 
V <19 (M/MQ)'''^*km s~^. Nonetheless, in view of the large observational uncertainties, along with the fact that the bulk of 
the ADAF emission originates from < lOOrg, the models should also be suitable for even higher velocity black holes. 
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With these results in hand, we may assess the number of RBHs that are potentially detectable by 
near future instruments. As the nearest example of a high density region of the ISM, we take the Orion 
molecular cloud, with a size of ~ 50 pc, density ngas lO'^ cm~'^ and distance ^ 400 pc (Lada, Strom, & 
Myers 1993). RBHs in this cloud would generate X-ray luminosities in excess of ^ 10^^ erg s~^ at ^ 4 keV 
if they move at <^ 20 km (Figure 1); this translates into fluxes ~ 5 x 10~^^ erg s~^ cm~^, which 
are above or near the expected detection limits of AXAF, XMM and ASTRO-E. Assuming RBHs obey a 
Maxwellian velocity distribution, the fraction of RBHs moving slower than V is 



The present velocity dispersion a of RBHs is an unknown quantity, and depends on the their uncertain 
initial velocities and positions (Campana & Pardi 1993). Here we simply take cr = 30 km s^^, the velocity 
dispersion of population I stars, which results in /(< 20 km s^^) = 0.07. The average number density of 
IOMq RBHs in the solar neighborhood has been estimated in a simplified manner by Shapiro & Teukolsky 
(1983) to be 8 X 10"'' pc~'^. With this value, the Orion molecular cloud would harbor a total of ~100 
RBHs, of which 7 could be detected by the aforementioned X-ray observatories. More detailed studies 
of RBH formation and evolution have been conducted by Campana & Pardi (1993), and if their somewhat 
model-dependent results are used, the expected number becomes <^ 6. 

The most slowly moving RBHs may also be observable even if they do not reside in molecular clouds. 
RBHs with V = A km s^^ in the average ISM density of rigas = 1 cm"*^ can radiate as much as those with 
V — 20 km s~^ and rigas = 10^ cm^ (Eq. (2|). Thus the same X-ray instruments as above are capable 
of detecting RBHs in the general ISM within 400pc of the Sun. With /(< 4 km s"') = 6.3 x lO""* for 
cr = 30 km s~' and the thickness of Galactic disk being 100 pc, we obtain in total ^ 25 such RBHs. 

Observations at lower frequencies should also play a role in the identification of RBHs, as their power 
output reaches a maximum in the infrared-optical region. The expected mid-IR flux from RBHs at the 
distance of Orion arc 0.1 mJy at ~ 10 fim for V <^ 50 km s^', above the typical detection thresholds of 
telescopes such as IRIS or SIRTF. Wide-field surveys by these facilities planned to hunt for brown dwarfs 
could uncover a sizable number of RBHs as a byproduct. Heavy extinction may completely obscure the 
optical to UV emission from RBHs within molecular clouds, and should reduce their near-IR and soft X-ray 
fluxes as well, although these wavebands may still be accessible for RBHs outside dense parts of the ISM. 

To discriminate RBHs from other IR or optical sources (e.g. young stellar objects), one calls for more 
information than narrow band spectra. Several effects could be of help: variability on various timescales 
from <^ 1 day to ~ 1 yr, possibly arising from fluctuations in accretion rate caused by ISM inhomogeneities; 
locational distribution, expected to be more or less uniform throughout our neighborhood, as opposed to 
young stars preferentially inhabiting molecular cloud cores with small filling factors (< 0.01); and proper 
motions which may be detectable for the more closer ones. Yet the most powerful method of identifying 
RBHs should be combining multiwavelength data. A preferred strategy may be to search first for candidate 
hard X-ray sources, where there is the least worry over interstellar absorption or confusing objects, and 
then conduct radio, IR and optical follow-up observations to see if the anticipated spectral distribution is 
revealed; this should prove effective in separating RBHs from other X-ray sources. 

We point out here one theoretical uncertainty regarding the low frequency emission from ADAFs. Our 
models assume the particles are always thermal, but Mahadevan & Quataert (1997) have recently shown 
that ADAFs at accretion rates M <^ 0.4a^MEdd contain regions where the electrons do not thermalize, their 
energy distribution being determined by the relevant heating and cooling processes. At AI <^ lO^'^a^MEdd 




(3) 



- 5 - 



this applies throughout the flow, in which case a nonthermal distribution with a reduced fraction of particles 
in the high-energy tail may be realized, leading to large decreases in the synchrotron and Compton emission 
components. The consequences should not be so pronounced, however, for emission components that are 
caused mainly by electrons with the plasma's mean energy, such as bremsstrahlung (indicated separately for 
one case in Figure 1); thus the hard X-ray fluxes, critical to the present discussion, are not greatly modified 
by this effect. 

ADAFs may also emit high energy gamma-rays at luminosities comparable to the bremsstrahlung 
component as a result of pion production in collisions between hot protons (Mahadevan, Narayan & Krolik 
1997). Despite the known association of some unidentified EGRET sources with high density regions of the 
ISM (e.g. Yadigaroglu & Romani 1997), the gamma-ray luminosities expected from RBHs are probably 
too low to account for any of these sources. Nevertheless, next generation instruments like GLAST may be 
sensitive enough to detect gamma-ray emission from the nearest RBHs. 

The above conclusions are in marked contrast to previous studies (McDowell 1985, Campana & 
Pardi 1993, Heckler & Kolb 1996) with the IP model, the results of which are overlayed in Figure 1 for 
comparison. The IP spherical accretion model is based on the ad hoc assumption of equipartition between 
gravitational, kinetic and magnetic energies, and possesses a single, high temperature, which results in 
synchrotron emission chiefly in the IR to soft X-ray range. At the pertinent accretion rates the hard X-ray 
flux from Comptonization in this model (Ipser & Price 1983) is insigniflcant, and as UV to soft X-rays suffer 
severe interstellar attenuation, visibility is practically limited to the IR-optical range. As noted above, 
IR-optical spectral data is insufficient for unambiguously distinguishing RBHs from other sources in the 
same waveband which are far more numerous, and one must rely on difficult measurements of variability 
or proper motion. RBHs accreting via ADAFs should be much easier to distinguish by their high energy 
emission. 

Heckler & Kolb (1996) have recently suggested that the color-based source selection scheme and spectral 
measurements intended for QSO searches in the Sloan Digital Sky Survey (SDSS) may also be advantageous 
for discerning isolated black holes that radiate a la Ipser & Price. Although ADAFs predict optical colors 
notably different from the IP model, and fluxes which are lower by about two orders of magnitude (and 
perhaps even lower in the case of a nonthermal electron distribution). Figure 1 suggests that the X-ray 
detectable RBHs are also visible to the SDSS, if extinction is not serious. With its exceptionally large fleld 
of view, the SDSS could still be useful if its results are utilized in conjunction with X-ray information. 

3. Emission from MACHOBHs 

As remarked in the introduction, there is the intriguing possibility that a large number of low-mass 
black holes comprises the dark halo of our Galaxy and causes the microlensing events. Here we choose 
a MACHOBH mass of 0.5Mq, the average value implied by microlensing event durations. The emergent 
spectra in the ADAF model have been calculated for such MACHOBHs with the same procedure as for 
RBHs and are shown in Figure 2. 

The smaller mass and lower accretion rate (Eq. ||^) compared to RBHs make the overall luminosities 
of MACHOBHs substantially lower. For a medium of rigas = 10^ cm^"^, 400 pc away, observability demands 
velocities V <^20 km s~^ for the near- IR band, and V ^ 4 km s^^ for X-rays. The fraction of MACHOBHs 
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whose velocities are smaller than V is given by 



V /.27r 



/ma(< V^) = xo/a 3 / v'^dv / dcj) d6'sin^exp( — ^ 



+ Kot + 2?;V^rot cose* 



MA -'O ^0 JO \ ^"MA 



exp ( --f^ / dt;^;sinh exp --^ • (4) 

TT Krot^MA V yiK J h K'^MaJ \ ^CTmA/ 

Here ctma is the MACHO velocity dispersion, and Vrot is the Galactic rotation velocity in the solar 
neighborhood with the direction of rotation being {9,(j)) = (0,0). We take ctma — 155 km s^^ 
and Viot ~ 220 km s^^, and assume the density of dark matter in the solar neighborhood to be 
0.0079 Mq pc-3 (Alcock et al. 1997), which gives a MACHOBH number density of 0.016 pc'^. As 
/ma(< 20 km s~^) = 2 X 10~^, the number of near-IR observable objects in Orion is only 0.4, while the 
chances of an X-ray detection are hopelessly small. For MACHOBHs in the general ISM of Ugas = 1 cm~^ 
the detection requirements are V < 4 km s"-'^ (near-IR) and V < 0.9 km s^^ (X-ray) at distances ^ 400 
pc; /ma(< 4 km s^^) — 2 x 10^® then results in ~ 2 near-IR MACHOBHs, with the X-ray numbers again 
being miniscule. 

MACHOBH spectra in the IP spherical accretion model are also shown in Figure 2. Were the model 
applicable, its greater luminosity at optical wavelengths may raise the number of detections to ~ 10. 
As remarked in the previous section, however, discriminating isolated, accreting black holes from other 
sources based solely on IR-optical observations is akin to searching for needles in a haystack. Applying 
either the ADAF model or the IP model, we must conclude that identification of MACHOBHs by their 
accretion-driven emission is extremely difficult if not impossible. 



4. Conclusions 

Implementing the advection dominated accretion flow model, we have considered accretion from the 
interstellar medium and the resultant emission for black holes of two types, those originating as stellar 
remnants and those constituting MACHOs. A fair number of stellar remnant black holes may be observed 
by X-ray, infrared and optical telescopes in the near future, which would provide strong proof of their 
existence. MACHO black holes, due to their smaller mass and larger velocity on average, are probably too 
underluminous for clear identification; gravitational radiation from binary MACHOBH mergers may be the 
only means to see them directly in the foreseeable future (Nakamura et al. 1997b). 

We may ask one final question: could accretion-induced radiation be of any significance for the case of 
neutron star (NS) or white dwarf (WD) MACHOs? Unlike black holes, radiation from accretion onto a NS 
or a WD is dominated by thermal emission from their surfaces where the radiative efficiency is high, and 
ambiguities in the model of accretion is of no serious concern. Such objects accreting from the ISM should 
indeed be capable of producing observable radiation if they are slow enough. The main difficulty we face 
here is misidentification with objects evolved from disk stars: simple estimates demonstrate that NS or WD 
MACHOs which meet the criterion for detectability are outnumbered by objects of disk origin with similar 
properties. Thus identifying them by their accretion-induced emission is also problematic. 

We thank H. Hanami, F. Iwamuro, Y. Kan-ya, S. Mineshige, T. Ohashi, S. Oya, M. Sakano and N. 
Yamasaki for helpful discussions. We are also grateful to an anonymous referee for several suggestions 
which improved the paper. This work was supported in part by the JSPS Research Fellowship for Young 
Scientists, and Grant-in- Aid of the Ministry of Education, Culture, and Sports No. 09640351. 
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Fig. 1. — The luminosity spectra of RBHs in the ADAF model (solid hnes) with M = 10 Mq and four 

different values of M (marked on each curve in Eddington units), which correspond to black hole velocities 
V = 10,20,40,80 km s^^ for rigas = 10^ cm^^^ or F = 2.2,4.3,8.6,17.2 km s"! for ngas = 1 cm'^. The 
bremsstrahlung component alone is exemplified by the dotted line in one case (M = 6.7 x 10~^MEdd)- To 
be compared are the Ipser-Price model spectra (dashes lines) for two M values. The typical point source 
detection limits for various observational facilities mentioned in the text have also been indicated, assuming 
a distance of 400pc (e.g. Orion). 
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Fig. 2.— The luminosity spectra of MACHOBHs in the ADAF model (solid lines) with M = 0.5 Mq and 
three different values of M, corresponding to ^ = 5, 10, 20 km s~^ for rigas = 10^ cm~^. Also shown are the 
Ipser-Price model spectra (dashes lines) for two M values. Otherwise similar to Fig.l. 



